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ABSTRACT 

We analyze the phase diagram of compact QED on the torus with a chirally 
symmetric four fermion interaction. Inside a mean field approximation for 
the four fermion term, a line of first order phase transitions and another one 
of second order are found in the {P,G) plane. Approaching the second order 
line a continuum limit can be defined. Critical exponents vary along this line 
in a similar way as in the non-compact model, suggesting that a non trivial 
interacting continuum theory can be constructed. 
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1. Introduction 

Non perturbative analysis of strongly coupled QED started several years 
ago. The aim of these investigations was and is the search for non gaussian 
fixed points in the abelian gauge model; the main motivation for this kind of 
business being the relevance of non asymptotically free gauge theories in the 
context of composite models for the dynamical mass generation mechanism 
in the Standard Model. 

First attempts to understand the nature of the continuum limit in full 
QED with dynamical fermions were done using the standard Wilson com- 
pact action for the pure gauge term. However there is no reason to prefer the 
compact formulation against the non compact one which, as well known, pre- 
serves also gauge invariance in the lattice formulation of the abelian model. 
Furthermore it was also argued that the non compact formulation, being free 
of monopoles as local minima of the action, would be nearest to the con- 
tinuum formulation. This beside the early finding that the compact model 
suffers from a strong first order phase transition in the chiral limit ad- 
dressed the attention to the non compact formulation of the model where 
evidence for the existence of a non gaussian fixed point has been found |Q . 

There are also several more recent attempts to reexamine the old re- 
sults on the first order character 0, § of the phase transition in the pure 
gauge compact model. All these attempts are based in the suppression of 
monopoles well by modifying the standard Wilson action with the inclusion 
of a monopole term [^] or by changing the standard torus topology [0, |]. 
In both approaches, which can be regarded as interpolations between the 
compact and non compact formulations but were the underlying pure gauge 
theory is not quadratic in the gauge field, the results suggest the existence of 
a non gaussian fixed point. Unfortunately no results with dynamical fermions 
are yet available here. 

We want to show in this paper how the compact formulation of the abelian 
model with a chirally symmetric four fermion interaction has a line of con- 
tinuous phase transition points, where a continuum limit can be defined. We 
will also show in the context of a mean field approximation for the fermion 
field that critical exponents take non mean field values, at least for critical 
points corresponding to not very small values of the bare gauge coupling. 

The introduction of the four fermion term in the action is not an ar- 
bitrary attempt to modify the known results of the compact model. In 
fact there are strong arguments based on a non perturbative analysis of the 
Schwinger-Dyson equations in the continuum formulation |Q suggesting that 
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four fermion interactions become renormalizable outside perturbation the- 
ory. This is not surprising at all since naive dimensional arguments, which 
would exclude four fermion interactions as renormalizable, should not work 
near non gaussian fixed points where large anomalous dimensions can be 
generated. 

The same analysis we are going to present here has also been done for 



the non compact formulation of the model [10 . As we will show, no quali 



tative differences in the phase diagram of both models appear in the weak 
gauge coupling regime, a region were we expect monopoles are dynamically 
suppressed. 

The gauged Nambu-Jona Lasinio (GNJL) model regularized on a lattice 
and with compact gauge fields has the following action. 



S = -P Yj cosQ^^{n) + ijA{9)tp + mi)^p - G ^ V„^„^„+^Vn+/, (1) 

where f3 is the inverse square gauge coupling, A(^) the massless Dirac oper- 
ator for Kogut-Susskind fermions, G the four fermion coupling and 

Q^n) = 9^{n) + e,{n + fi) - e^{n + 0) - 9,{n) (2) 

In the chiral limit (m = 0) the action (1) is invariant under the continuous 
transformations 



which define a.U{l) symmetry group. 

Since the four fermion interaction is not bilinear in the fermion field, it is 
not suitable for numerical simulations. To overcome this problem, we need 
to introduce an auxiliary vector field interacting with fermions, the path in- 
tegral of which reproduces the four fermion coupling. To perform numerical 
simulations of such a system is highly complex since it contains three pa- 
rameters and G) and two vector fields [|1T|. Alternatively, we can use a 



mean field approximation to bilinearize the fermion action. Although quan- 
titative results will be affected by this approximation, we expect qualitative 
conclusions to be realistic . 



As usual in the standard mean field technique, we make in (1) the follow- 
ing substitution: 
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Gj2i'ni^ni'n+f,i'n+fi 2dG (4)^: ) ^ i^n'ipn (4) 

n,ii n 

where d is the space-time dimension. The action (1) becomes in this way a 
bihnear in the fermion field and the path integral over the fermion degrees of 
freedom can be performed by means of the Matthews-Salam formula. In the 
following sections we shall analyze the phase structure and critical behavior 
of the model within this mean field approximation. 



2. The Phase diagram 



The analysis of the phase diagram is similar to that of ref. |I0[, and 
details can be found there. We shall study the chiral condensate and the 
plaquette energy. The chiral condensate is an order parameter for chiral 
symmetry breaking. Within the mean field approximation (4) it is given by 
the selfconsistent equation: 

where V is the lattice volume and ±«Aj (j = 1, . . . , \^/2) are the eigenvalues 
of A, which depend on the compact gauge field 9. The expectation value 
of the r.h.s. of (5) is taken over the compact gauge field configurations and 
includes in the probability measure, besides the exponential of the pure gauge 
action, the fermionic determinant of standard QED evaluated at the effective 
mass ffi = m — 8G{ ipip ) . 

In the chiral limit, eq. (5) becomes 

(,.,>.16G(^^)(i|^^^^jI— ,) (6) 

This equation is always verified if (ipip) = 0, and this is the only solution 
in the symmetric phase. In the broken phase, where (ipip) ^ 0, the chiral 
condensate is given by the equation: 

l = 16G'(-!-V _ , ) (7) 
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the solution of which minimizes the free energy. It is therefore the physical 
(stable) solution. In the G — > oo limit the solution of eq. (7) at leading order 
is 

1 = 16G ^—^ + 0(-^) (8) 

This shows that at large G chiral symmetry is spontaneously broken for all 
values of /5, the order parameter being 

(^^)~^ (9) 

The existence of a second order phase transition line separating the bro- 
ken and symmetric phases follows from eq. (7) |]10[. It gives the following 
expression for the critical line: 

^^(/^) = ia,J/2 1. (10) 

where the probability distribution function in the previous mean value is 
that of massless compact QED. The r.h.s. of eq. (10) is proportional to the 
inverse transverse chiral susceptibility of compact QED in the chiral limit, 
which can be obtained from the identity 

This allows us to discuss qualitatively the phase diagram of the compact 
GNJL model from the known properties of compact QED. 

First, let us remember the main features of compact QED. In the chiral 
limit a strong first order phase transition takes place at a value of the coupling 



/3c ~ 0.9 ||T2|- It separates a strongly coupled phase, (5 < Pc, where chiral 
symmetry is spontaneously broken, from a Coulomb phase, /? > Pc, char- 
acterized by a symmetric vacuum. At (3 = f3c the chiral condensate shows 
a jump from zero to some non-zero value, and also the plaquette energy 
jumps. This first order phase transition continues at non-zero bare fermion 
mass 1^, 13| along a line in the (/?, m) plane which connects the {Pc, 0) point 
in the chiral limit with the (/3q, oo) transition point of the quenched theory, 
at m = oo ( /9„ ~ 1 , p ) . Both the chiral condensate and plaquette energy 



are discontinuous along this line |T2|. We denote this line by (/5,mc(/?) 



Now, let us switch on the four fermion interaction. For P values larger 
than Pc the transverse susceptibility of compact QED is finite in the chiral 
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limit, and we get a line of second order phase transitions eq. (10) which 
separates a phase where chiral symmetry is spontaneously broken, when G > 
Gc{P), from a symmetric phase, when G < GdP). When /? approaches (3c 
from the right hand side, the QED massless transverse susceptibility remains 
finite, and thus GdPc) > 0. This is a significant difference respect to the non 
compact model where, due to the fact that the susceptibility diverges in the 
P ^ P+ limit, GciPc) = 0. 

In the broken phase of compact QED the transverse susceptibility is di- 
vergent in the chiral limit due to the Goldstone theorem, eq. (11). Therefore 
chiral symmetry is spontaneously broken for all values of G if /? < f3c- On 
the vertical line (3 = j3c, < G < GciPc) the transition is first order since in 
this range of G the chiral condensate jumps from a non vanishing value at 
/3 < /?c to zero at /5 > Pc- The plaquette energy is also discontinuous along 
this line. 

The first order line at m 7^ in compact QED results in the chiral limit 
of the compact GNJL model in a first order phase transition line which 
prolongates the vertical line from the point {Pc,Gc{,Pc) ), ending at Pq, G = 
00. These first order transitions are originated by the effective fermion mass 
m generated when ( ipip ) 7^ 0. Since the plaquette energy and the chiral 
condensate for the compact GNJL model in this mean field approximation 
are obtained from those of compact QED evaluated at an effective fermion 
mass m, it is clear that for all Pc < P < Pq there exists a G/o > GdPc) such 
that 

-Gfo{^^){P,Gfo)=mc{P) (12) 

with Gfo(Pg) = 00. On this line {P,Gfo{P)) the plaquette energy and the 
chiral condensate show a jump, producing a first order phase transition. In 
Figure 1 we show a qualitative picture of the phase diagram for the compact 
GNJL model. 

3. Critical exponents 

Since second order phase transitions are usually related to divergent cor- 
relation lengths, we expect that a continuum limit can be defined along the 
critical line given by eq. (10). The nature of the renormalized theory will 
depend on the critical exponents at each critical point. The calculation of the 
critical exponents in this compact model is analogous to that of the noncom- 
pact version [|l^]. Here we will briefly comment the way to calculate them 
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and refer the reader to ref. |T0| for details. 

We write the chiral condensate of compact QED (without the four fermion 
couphng) as 

= -2mF{p,m) (13) 
where F is a function of jS and m, whose behavior when m ~ is 

F{(3,m) = F{(3,0) + Em" + ... (14) 

Equation (13) evaluated at the effective fermion mass m = m—'^G^ipip ) is 
the equation of state for the GNJL model in our mean field approach. From 



(13) and (14) and after simple algebra |]T0| we get the following relations for 
the critical exponents: 

5 = u + l Pm = - 7 = 1 (15) 

UJ 

The practical rule to measure the critical exponents in the GNJL model 
at a given critical point (/5, G(/3)) is therefore to measure the u exponent of 
the subleading contribution to the chiral condensate in the Coulomb phase 
of compact QED. This exponent can be related to the behavior at the origin 
of the spectral density of eigenvalues of the Dirac operator. In fact, writing 

— /'°° 2,171 

(V'V^) = -/ d\p{\)——^ (16) 
the following relation is derived: 

uj = p-l (l<p<3) to = 2 ip>3) (17) 

where p in (17) is the exponent which controls the behaviour of p(A) near the 
origin, i.e., p(A) ~ A^ when A — 0. Notice that p is always larger than 1 in 
the Coulomb phase of QED as follows from a finite massless susceptibility. 

The value of p can be obtained from the analysis of what we call gener- 
alized susceptibilities, defined by 




X. = T7 Et^ • (18) 



They are finite if g < p + 1 and divergent otherwise. On finite lattices we 
expect 

X, ~ A, + 5,L-('^-P-i) (19) 
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with a logarithmic divergency with the lattice size in the especial case q = 
p+1. Ag and Bg in (19) are constants independent of the lattice size and a is 
some positive number. These finite size scaling relations allow us to calculate 
p and therefore the critical exponents. 

The procedure we use is as follows. From simulations of compact QED 
in the chiral limit using the MFA approach |0, we are able to compute the 
generalized susceptibihties for several values of q and lattice size L. Fitting 
them with the function 



X, 



ag + hnL~ 



(20) 



it is possible to extrapolate the results to the thermodynamical limit. We 
obtain in this way the critical value qc such that > for q < qc and ag = 
for q > qc- Eq. (19) implies p = qc — I, and from (17) we get 
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P. 



(21) 



gc-2 

We have applied this procedure to two points on the critical line of Fig. 1, 
(3 = 1.020, G = 0.078 and (3 = 1.045, G = 0.099. The corresponding critical 
values of q have been obtained from the finite size scaling analysis (20) on 
4'*, 6"^, 8"^ and 10^ lattices. The extracted qc values together the corresponding 
critical exponents are reported in Table I. As in the non compact GNJL 
model [10 1, critical indices far from their mean field values are obtained for 



(3 values near the phase transition point of QED. 





Gc 


Qc 


Pm 


S 


7 


1.020 


0.078(1) 


1.5(2) 


2.0(6) 


1.5(2) 


1 


1.045 


0.099(3) 


2.7(2) 


0.59(7) 


2.7(2) 


1 



4. Discussion 



We have derived in this paper the phase diagram of the compact GNJL 
model. A first order line starting from the critical coupling of the QED 
model appears for every value of the four fermion coupling, ending at the 
critical value of quenched QED at G = oo. What is more exciting however 
is the finding of a second order line in the compact formulation of the model. 
This continuous phase transition (discontinuous line in Fig. 1) prolongates to 
the weak gauge coupling region approaching the P = oo axis at some finite 
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value of G, a region where certainly monopole condensation or monopole 
percolation does not hold. 

We have also shown how critical exponents, computed in the mean field 
approximation for the fermion field but taking into account fiuctuations of 
the gauge field, take non mean field values and change along the critical line. 
This is very similar to what we found in the non compact version of the 



model [|TOl- Again in the compact formulation, fiuctuations of the gauge field 
seem to play a fundamental role in driving critical exponents to non gaussian 
values. One is then tempted to speculate that both formulations would drive 
to the same renormalized continuum physics. 

The last but not less important feature of the compact formulation of 
the gauged Nambu-Jona Lasinio model is the relevance of the four-fermion 
coupling in the search for a continuum limit. In contrast to what happens 
in the non compact formulation the four-fermion coupling, that as stated 
before becomes renormalizable outside perturbation theory, is an essential 
ingredient to define a continuum limit for the compact model. 
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Figure captions 



Figure 1. Phase diagram of the compact GNJL model in the {P,G) plane. 
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